Electronic signature of DNA nucleotides via transverse transport 
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We report theoretical studies of charge transport in single-stranded DNA in the direction per- 
pendicular to the backbone axis. We find that, if the electrodes which sandwich the DNA have the 
appropriate spatial width, each nucleotide carries a unique signature due to the different electronic 
and chemical structure of the four bases. This signature is independent of the nearest-neighbor nu- 
cleotides. Furthermore, except for the nucleotides with Guanine and Cytosine bases, we find that the 
difference in conductance of the nucleotides is large for most orientations of the bases with respect 
to the electrodes. By exploiting these differences it may be possible to sequence single-stranded 
DNA by scanning its length with conducting probes. 

PACS numbers: 73.63.-b, 87.15.-v, 87.64.-t 



There has recently been an increased interest in charge 
transport in DNA, due to both its relevance in physi- 
ological reactions and to its potential use in molecular 
electronics. 0, 0, HJ Previous studies have looked into 
the effect of the base sequence and structural distortions 
on charge transport, and the interplay between different 
transport mechanisms. Q H S HE Q, H U ES El 03 
However, much of the research so far has focused on how 
charge flows along the DNA axis. Very few experimental 
studies have looked into the transport properties of DNA 
in the transverse direction. For instance, scanning probe 
techniques have been used to image isolated DNA bases 
or a few bases in single-stranded DNA [HE1EEEEE3 
but no systematic study is available on the relative ability 
of each base to carry electrical current. 

In this paper, we report on charge transport in single- 
stranded DNA in the direction perpendicular to the back- 
bone axis. Our goal is to determine the relative cur- 
rent carried by each nucleotide in the presence of its 
neighbors. We envision a possible experimental setup 
as schematically shown in Fig. ^ where two probes of 
nanoscale dimensions sandwich the DNA in such a way 
that one or a few bases are spatially confined within the 
lateral dimension of the probes. Experimentally, this 
could be realized by a polynucleotide molecule between 
two electrodes embedded in a nanopore 0, E3 or by ad- 
sorbing the polynucleotide onto a metal surface and using 
a scanning tunneling microscope (ST M) tip t o measure 
the current through a single nucleotide. |14ill5| We there- 
fore seek to identify specific electronic signatures of each 
nucleotide that would allow us to create a DNA elec- 
tronic map. In particular, we investigate the role the 
electronic and chemical structures of the bases play in 
the charge transport, as well as the role of interactions 
between each base and its nearest neighbors. We indeed 
find that the current through each individual nucleotide 
provides a unique electrical signature of each base. The 



ratio of the current for different nucleotides gives a mea- 
sure of the difference of their signatures. These electri- 
cal signatures change with orientation of the nucleotides, 
but are only affected slightly by neighboring nucleotides 
so long as the electrode lateral dimension is comparable 
with the average "size" of the bases (i.e., a dimension on 
the order of 1.0 nm, comparable to the nucleotide spac- 
ing of 0.7-0.9 nm [ljj). In most cases, the current ra- 
tios are quite large when comparing between nucleotides 
of random orientations, except for the ratio of the cur- 
rent through 2'-deoxyguanosine 5'-monophosphate and 
the current through 2'-deoxycytidine 5'-monophosphate 
(i.e., nucleotides with guanine and cytosine bases). Such 
findings may help create a practical sequencing tool based 
on the electronic read-out of DNA. 

We consider a single strand of DNA, e.g., a polynu- 
cleotide, between two electrodes kept at a given bias, 
(see schematic in Figure . The electrodes are repre- 
sented by 3 layers of about 30 atoms arranged in the 
[111] surface geometry, which gives an electrode surface 
spanning about 100 A 2 . The electrode-electrode spacing 
is fixed at 15 A. Both the area of the electrode surface 
and the electrode-electrode distance are such that the 
largest of the nucleotides (i.e., the one with the A base) 
can be accommodated between the electrodes. Since we 
are interested only in the ratio between currents of dif- 
ferent nucleotides (and we consider only linear response) , 
we employ a simple tight-binding model for the elec- 
tronic structure of the system with s,p x ,p y , and p z or- 
bitals for each carbon, nitrogen, oxygen, and phosphorus 
atom; and s orbitals for hydrogen and gold. We then 
use a Green's function approach in order to calculate 
the current through the whole system (electrodes plus 
DNA). H^IIJ The Green's function Q DNA of the com- 
bined electrode-DNA-electrode system is written as 
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FIG. 1: Schematic of a polynucleotide between two electrodes. 
The electrodes only couple to a single nucleotide at a time. 



where Sdna and TLdna are the overlap and the Hamil- 
tonian matrices, respectively. ^t(b) are the self-energy 
terms that describe the effect of the electronic structure 
of the leads. 

The transmission coefficient T(E) is given by [2(J 

t(e) = T a = Tr\r t g DffA r b gl, NA ] (2) 

where ^t(b) = — We calculate the current 

from 

I= T j dET(E)[f t (E)-f b (E)}, (3) 

where f t (b) = {exp[(£' - fi{(ij)/feT] + is the Fermi 
function, and fJ-ttb) 1S the chemical potential of the top 
(bottom) electrode (see Fig^l. We assume the voltage 
drops equally in the space between the top electrode and 
the molecule, and the molecule and bottom electrode. 
The current is calculated at room temperature. 

We first examine the current ratios for the four nu- 
cleotides (which we label by A, G, C, or T corresponding 
to their base) without nearest-neighbor nucleotides at- 
tached. The nucleotides are therefore "passivated" , i.e., 
the oxygen atoms which are normally charged in solu- 
tion are instead bonded to a hydrogen atom. [l| For each 
nucleotide, we relaxed the geometry using total-energy 
Hartree-Fock calculations. |23j We consider first the nu- 
cleotides positioned so that the base is in a plane per- 
pendicular to the electrode surface (see Fig. 0). Figure |2 
shows the ratio between the current of A and the current 
of the other nucleotides, Ia/Ix where X = G, C, or T. 
The ratios are about 20, 40, and 660 for G, C, and T, 
respectively. Strong contrast has indeed been observed 



in an STM experiment between A and T nucleotides. [l4j 
Figure 2 thus indicates that the isolated nucleotides can 
be distinguished relatively easily from one another using 
transverse currents. 

The differences in the electrical current ratios are due 
to the relative position of the Fermi level with respect to 
the highest-occupied molecular orbital (HOMO) and the 
lowest-unoccupied molecular orbital (LUMO) as well as 
the value of the density of states (DOS) |24j of the base at 
the Fermi level. By projecting onto the different atomic 
orbitals, we find that the states of the bases make up 
for most of the HOMO and LUMO levels while the DOS 
at the Fermi level is mainly determined by the backbone 
states due to the large coupling of the backbone with the 
bottom electrode. However, the relative current is closely 
related to the value of the DOS at the Fermi level due 
to the bases, which gives the ordering A > G ~ C > 
T. How well those base states couple to both electrodes 
determines the overall magnitude of the relative currents. 

Having shown that each nucleotide carries a different 
electrical signature for a particular geometry, we can now 
examine how this result changes with orientation of the 
base with respect to the probe surfaces and location of 
the backbone atoms relative to the surface. We define 
our coordinate system as shown in Figure|21 The origin is 
the nitrogen which attaches each base to the phosphate- 
sugar backbone. Both rotation and translation involve a 
rigid movement of the whole nucleotide (base plus back- 
bone). The x-axis runs parallel to the electrode surfaces 
and perpendicular to the plane of the base. The y-axis 
runs parallel to the electrode surfaces and in the plane of 
the base. The z-axis runs perpendicular to the electrode 
surface. We consider rotations and translations of the 
DNA nucleotides about each of the axes independently. 
The six variations we consider are the following: (1) ro- 
tation about the x-axis from -10 to +10 degrees, (2) rota- 
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FIG. 3: Coordinate system which we use to define rotations 
and translations of the nucleotides with respect to the probe 
surfaces. The x-axis is pointing out of the paper. The dashed 
lines represent the position of the electrode surfaces. 




FIG. 4: Current ratios of A to the other nucleotides X (X — 
G, C, or T) at a bias of 0.1 V . I a is for A at its original config- 
uration, while the other nucleotides' orientations are varied. 
The 6 lines for each nucleotide, from left to right, correspond 
to rotation about the x-, y-, and z-axes and translation in the 
x-, y-, and z-directions. 



tion about the y-axis from -15 to 15 degrees, (3) rotation 
about the z-axis from -30 to +30 degrees, (4) translations 
along the x-axis by -1.5 to 1.5 A, (5) translations along 
the y-axis by -1.5 to 1.5 A, and (6) translation along 
the z-axis by -0.5 to 0.5 A about the original nucleotide 
position. |25| 

Figure 01 shows relative changes in current with nu- 
cleotide rotation and translation as defined above. For 
each nucleotide, the lines from left to right correspond to 



rotation about the x-, y-, and z-axes and to translation 
along the x-, y-, and z-axes, respectively. The variations 
(2) and (6) (i.e., rotation about the y-axis and trans- 
lation along the z-axis, respectively) cause the largest 
changes in relative current. This is easy to understand 
since both movements change the distance between the 
uppermost atoms on the nucleotides and the electrode 
surface, which gives a large change in coupling. This is 
analogous to the double-stranded DNA case: variation 
of the angle and distance between bases can significantly 
change the conductivity because of a reduction of the 
inter-base coupling. [12( From Fig. 0] it is clear that, 
even taking into account the above variations, there is 
a large distinction between the nucleotides A and T and 
the pair {G,C}. The latter two have a similar electrical 
signature for a wide range of rotations and translations 
and cannot therefore be easily distinguished. Shot noise, 
however, being more sensitive to changes in electronic 
structure than the current, j2|| could be used to distin- 
guish between the latter two bases. We leave this study 
for future investigations. 

We conclude by looking at the effect of nearest- 
neighbor nucleotides on the electrical signature of each 
individual nucleotide. It is known from ground-state cal- 
culations [IlillSIIil that in isolated stacks of nucleotides 
the states on each base are highly localized on that base. 
We therefore expect small corrections due to nearest- 
neighbor base interactions. We have checked this point 
by looking at the following representative sequences of 
three nucleotides of the type AAA, A ATT, CAA, GAG, 
GAT, and TAT, where A = G, C, or T is the nucleotide 
that is sandwiched between the two electrodes. |30j The 
base of the central nucleotide is also assumed perpendic- 
ular to the electrode surface. 

The relative change in current for these cases is shown 
in Figure |5] at 0.1V. As expected, there is little change in 
current compared to the individual nucleotides so long as 
the spatial width of the electrodes is only large enough 
to accommodate just one base (see Fig.^l. Most of the 
change is due to the difference in backbone geometry 
of the strand compared to the isolated nucleotide. Al- 
though, for sequences with T as the central nucleotide, 
the current changes due to nearest neighbors: T carries 
the smallest current and therefore is the most sensitive 
to its nearest neighbors interacting with the electrodes. 
We therefore conclude that the electrode width needs to 
be about 1 nm if the nearest neighbor bases have un- 
corrected orientations. When this is the case, there are 
variations due to single nucleotide orientation, as shown 
in Figure 01 but there is not significant variations due 
nearest neighbor interaction. 

In conclusion, we have studied charge transport 
through nucleotides in the direction transverse to the 
backbone axis. We find that each nucleotide has a 
unique electrical signature determined by the electronic 
and chemical properties of the bases irrespective of the 
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FIG. 5: Current ratios for strands of nucleotides with central 
nucleotide X (X — G, C, or T) at a bias of 0.1V. Ia is for 
isolated A at its original configuration. The crosses are for 
the individual nucleotides and the circles are for sequences of 
three nucleotides (see text for details). 



nearest-neighbor nucleotides. Except for the nucleotides 
with Guanine and Cytosine bases, these signatures are 
quite robust with variation in nucleotide orientation and 
location with respect to the electrodes. By exploiting 
these differences, it may be possible to sequence DNA 
by scanning its length with a conducting probe. Clearly, 
for such application many other factors (like solvent and 
ionic effects) need to be studied in detail and additional 
theoretical work is necessary. Also, the full noise spec- 
trum may provide extra diagnostic capability to differen- 
tiate between bases. 
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